Relaxation to non-equilibrium in expanding ultracold neutral plasmas 
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We investigate the strongly correlated ion dynamics and the degree of coupling achievable in the 
evolution of freely expanding ultracold neutral plasmas. We demonstrate that the ionic Coulomb 
coupling parameter ]?; increases considerably in later stages of the expansion, reaching the strongly 
coupled regime despite the well-known initial drop of I\ to order unity due to disorder-induced 
heating. Furthermore, we formulate a suitable measure of correlation and show that Ti calculated 
from the ionic temperature and density reflects the degree of order in the system if it is sufficiently 
close to a quasisteady state. At later times, however, the expansion of the plasma cloud becomes 
faster than the relaxation of correlations, and the system does not reach thermodynamic equilibrium 
anymore. 



PACS numbers: 52.27.Gr,32.80.Pj,05.70.Ln 

Freely expanding ultracold neutral plasmas (UNPs) 
1[ have attracted wide attention both experimentally 
1 H H H and theoretically H 0, H H 0. A main 
motivation of the early experiments was the creation of a 
strongly coupled plasma, with the Coulomb coupling pa- 
P-f rameter (CCP) V = e 2 /(ak B T) > 1 (where T is temper- 
ature and a is the Wigner-Seitz radius). From the exper- 
imental setup of Q , the CCPs of electrons and ions were 
estimated to be of the orders of r c w 30 and P; w 30000, 
respectively. By changing the frequency of the ionizing 
laser, the electronic temperature can be varied, offering 
the prospect of controlling the coupling strength of the 
electrons and creating UNPs where either one, namely 
the ionic, or both components could be strongly coupled. 

However, due to unavoidable heating effects |rj HH H^] 
these hopes have not materialized yet, and only r o « 0.2 
and T; « 2 have been confirmed. Furthermore, the evo- 
lution of the expanding plasma turns out to be a rather 
intricate problem of non-equilibrium plasma physics for 
which a clear definition of the degree of correlation is not 
obvious to begin with. 

The goal of this letter is twofold: Firstly, we will for- 
mulate a consistent measure of correlation for expand- 
ing ultracold plasmas, and secondly we demonstrate that 
the strongly correlated regime with Tj w 10 for the ionic 
plasma component can be reached by simply waiting un- 
til the plasma has (adiabatically) expanded long enough 
under already realized experimental conditions. This is 
remarkable in the l ight of alternatives proposed to in- 
crease Ti [l2l IT3L H4I flil Hil Il7| which are experimentally 
rather involved. 

Substantiating both of our statements theoretically re- 
quires the ability to propagate the plasma numerically 
over a long time with full account of the ionic correla- 
tions. To this end, we have developed a hybrid molecular 
dynamics (H-MD) method for the description of ultra- 
cold neutral plasmas. In our approach, ions and recom- 
bined atoms are propagated in the electronic mean-field 
potential with the full ion-ion interaction taken into ac- 



count. The much faster and weakly coupled electrons, on 
the other hand, are treated on a hydrodynamical level. 
Elastic as well as inelastic collisions, such as three-body 
recombination and electron-impact ionization, are incor- 
porated using a Monte-Carlo procedure 

[3E3- The 

H-MD approach accurately describes the strongly cou- 
pled ionic dynamics and therefore allows us to realisti- 
cally study the plasma relaxation behavior for long times. 

Assigning Ti for an expanding plasma by extracting a 
temperature from the kinetic energy of all ions is compli- 
cated by the fact that the radial expansion contributes 
considerably to this energy [lif. In our approach, we 
can determine a local temperature from the ion velocity 
components perpendicular to the (radial) plasma expan- 
sion . Additionally, the distribution of thermal veloci- 
ties of all plasma ions is found to be well described by a 
Maxwcll-Boltzmann distribution corresponding to an av- 
erage temperature Tj even at relatively early times. Ex- 
perimentally, the time evolution of the average ion tem- 
perature is determined from the correspon ding Do pplcr 
broadening of optical transition linewidths 19, 20]. The 
close agreement between experiment [2(j and theory (fig- 
ure^ supports both the experimental scheme of extract- 
ing an ionic temperature as well as the assignment of a 
temperature to the transversal ion velocities in the H-MD 
approach. 

Remarkably, the initial relaxation of the average ion 
temperature exhibits temporal oscillations, in contrast 
to the known behavior of weakly coupled plasmas. For 
the latter, the timescale t colT of the initial build-up of 
ion-ion correlations is typically much smaller than the 
timescale i rc i for the relaxation of the one-particle dis- 
tribution function. Based on this so-called Bogoliubov 
functional hypothesis, which is one of the fundamental 
concepts in kinetic theory [2]]], the different relaxation 
processes can be separated, resulting in a monotonic 
behavior of the correlation energy (and hence the ion 
temperature) j^. Molecular dynamics simulations of 
the relaxation behavior of homogeneous one-component 



FIG. 1: Calculated ion temperature (solid line) for a plasma 
of 10 6 Sr ions with initial peak density po(0) = 2-10 9 cm~ 3 and 
electron temperature T c (0) = 38K, compared to experimental 
results (dots) |2(I . The fact that the experimental ion number 
is about a factor of ten larger than in our calculation does 
not affect the time evolution of the ionic temperature, since 
there is no significant adiabatic ion cooling on the timescale 
considered in figure 



plasmas show that the ion temperature starts to un- 
dergo damped oscillations around its equilibrium value 
if both of these timescales become equal, which happens 
for Ti(0) > 0.5 (2^. Therefore, the nonmonotonic ion 
relaxation observed in ultracold plasmas may be seen as 
a direct manifestation of the violation of Bogoliubov's 
hypothesis. 

Compared to the homogeneous plasmas considered in 
[23| . the oscillations of the average ionic temperature 
damp out much quicker in the present case. This can 
be attributed to the fact that the Gaussian density pro- 
file of the UNPs created in current experiments leads to a 
spatial dependence of the correlation timescale i C orr, the 
build-up of correlations being fastest in the center of the 
plasma where the density is highest, and becoming slower 
towards the edge of the plasma cloud. As a consequence, 
the local ionic temperature shows not only temporal, but 
also pronounced spatial oscillations, which however tend 
to become averaged out if the spatial average over the 
whole plasma cloud is taken. 

Having established the approximate validity of assign- 
ing a global temperature to the plasma ions, it becomes 
possible to define a corresponding CCP 1Y While the 
initial ion relaxation reveals some interesting strong- 
coupling effects as discussed above, disorder-induced 
heating 0, Il2| drives the ion component to the border 
of the strongly coupled fluid regime I\ « 2 and therefore 
limits the amount of correlations achievable in UNPs. 
However, so far this could be verified only for the early 
stage of the plasma evolution |7], Il2l [l9j . The present 
H-MD approach allows us to study also the long-time 
behavior of the ion coupling. 

In figure |3 we show T; (solid line) as a function of 



20 40 60 80 100 

X 

FIG. 2: Ionic Coulomb coupling parameter for a plasma with 
M(0) = 5 ■ 10 4 , #(0) = 1.1 ■ 10 9 cm- 3 and T c (0) = 50K. The 
solid line shows the CCP calculated from the average temper- 
ature and density, the dashed line marks the CCP extracted 
from pair correlation functions (see text). Inset: blow-up of 
the short-time behavior. 

t = ojp.ot for a plasma with JVj(0) = 5 • 10 4 , Pi(0) = 
1.1 • 10 9 cm -3 and T c (0) = 50K, determined in a central 
sphere with a radius of twice the root-mean-square radius 
of the expanding plasma. (In the following, dimensionless 
units are used where time is scaled with the initial plasma 
frequency 07 p ,o = uj p (t — 0) and uj p = ^/4ire 2 pi/m,i.) 
As can be seen in the inset, T; quickly drops down to 
Ti i=a 2. After this initial stage, however, Ti starts to 
increase again due to the adiabatic cooling of the ions 
during the expansion. Indeed, CCPs of more than 10 are 
realized at later stages of the system evolution, showing 
that cold plasmas well within the strongly coupled regime 
are produced with the present type of experiments. 

Neglecting the influence of the changing correlation 
energy as well as inelastic processes, the adiabatic law 
for the plasma expansion yields T- 1 p i 2 ^ 3 = const.. 
Hence, Tj should increase ex p i 1 ' 3 as the plasma ex- 
pands, ultimately leading to coupling strengths of 10 2 or 
even larger at very long times. For a classical plasma in 
thermodynamical equilibrium, the Coulomb coupling pa- 
rameter is a direct measure of the amount of correlations, 
and properties such as pair correlation functions etc. can 
be parametrized by this single quantity. However, the 
UNPs created in the present type of experiments are non- 
equilibrium systems. Initially, e.g., they are created in a 
completely uncorrelated state, so that the high value of 
Ti caused by the ultralow temperature of the ions has no 
relation at all with the correlation properties of the sys- 
tem. At later times, the system relaxes towards a local 
equilibrium. However, the plasma is freely expanding, 
and hence constantly changing its steady state. Thus, 
the plasma is in a non-equilibrium state at all times, and 
one must ask to what extent Ti really parametrizes the 
correlations present in the plasma. 
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FIG. 3: "Pair correlation functions" of the plasma of fig. 
[2] at four different times r = 0.54, r = 1.1, r = 10 and 
t = 60.7. The Y\ indicated in the figure is obtained by fitting 
the distribution of scaled interionic distances (dots) with pair 
correlation functions for a homogeneous plasma given in |2q| 
(solid line). 



To this end, we compare T\ as obtained above with an 
alternative value Tj (dashed line in figure [2Jl parametriz- 
ing correlation properties of the plasma. As in |l5| . we 
have calculated the distribution P{r/a\ oc ) of interionic 
distances rescaled by the local Wigner radius. These dis- 
tribution functions are fitted to the known pair correla- 
tion function g{r/a,Y{) of an equilibrium plasma given 
in [2^ (figure 13 ■ From the fit, a value f; is extracted at 
several times. As can be seen in figure |3 at very early 
times the distribution of scaled interionic distances is not 
very well fitted to a pair correlation function of a homo- 
geneous plasma in equilibrium. Again, this is due to the 
fact that the system is far away from its steady state, 
and a single parameter does not describe the correlation 
properties of the plasma in an adequate way. However, 
the interionic distances quickly relax, and they are well 
described by a pair correlation function of an equilibrium 
system at later times. Hence, we conclude that the value 
of Tj is suitable for parametrizing the correlation prop- 
erties of the plasma cloud once it came sufficiently close 
to equilibrium, and that it indeed reflects the degree of 
coupling in the plasma. 

Comparing T; and Tj in figure [3 several conclusions 
can be drawn. As discussed above, and has been well 
known before, in the very early phase of the system evo- 
lution there is no relation between Tj and Tj since the 
plasma is too far away from equilibrium. As the plasma 
relaxes towards this equilibrium, Tj and T{ rapidly ap- 
proach each other, showing that during this stage Tj is a 
good measure for the correlation properties of the ions. 
In particular, the correlations building up in the system 
are indeed those of a strongly coupled plasma with a 
CCP well above unity. Moreover, the transient oscilla- 
tions characteristic of the relaxation process which are 



apparent in Tj also appear in Tj, however with a "phase 
shift" of 7r. This phase shift is due to the fact that a 
minimum in the temperature means a maximum in T; 
for a given density. Since total energy is conserved, a 
minimum in the thermal kinetic energy corresponds to 
a maximum in the potential energy, i.e. to an increased 
number of pairs of closely neighboring ions, and therefore 
to a pair correlation function with enhanced probability 
for small distances and consequently a minimum in Tj. 

At later times, both curves diverge again and the 
plasma evolves back towards an undercorrelated state. 
At first sight, this seems very surprising since the plasma 
should relax towards equilibrium rather than away from 
it. However, as argued above, the plasma is freely ex- 
panding and the corresponding equilibrium properties 
are constantly changing. We interpret figure [21 as being 
again evidence for the break-down of the Bogoliubov as- 
sumption of a separation of timescales, in this case of the 
correlation time r corr and the hydrodynamical timescale 
Thyd, i-e. the characteristic time for the plasma expansion. 

The timescale Th y d may be determined from the rela- 
tive change of macroscopic plasma parameters, such as 
the ion temperature or density. Due to the transient os- 
cillations of the ion temperature we choose the ion den- 
sity to characterize the change of the plasma properties 
(other choices such as, e.g., a cx p i 1 ' 3 lead to the same 
conclusions since they result in a simple constant propor- 
tionality factor 1 /a of order unity in the expression for 
T h yd)- Then 
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where we have used the selfsimilar solution for the colli- 
sionless quasineutral plasma expansion |2lj with r Gxp = 
cr(0)o;p ) o-\/rni/ (fcsT e ). On the other hand, binary correla- 
tions are known to relax on the timescale of the inverse of 
the plasma frequency in the strongly coupled regime |23j | 
for an initially uncorrelated state, and somewhat slower 
if the initial state already exhibits spatial ion correlations 
|l2j . Tcorr ^ ^p.o/^p- The selfsimilar plasma expansion 
then yields 



ft j 2/2 \3/4 
Tcorr = (1 + T /T c: I 



cxp/ 



(2) 



Therefore, r corr is initially much smaller than Th y d, but 
ultimately exceeds Th y d as the plasma expands, leading 
to an inevitable break-down of the Bogoliubov condition. 
Consequently, the build-up of correlations in the system 
cannot follow the changing equilibrium anymore, and cor- 
relations freeze out as indicated by the leveling-off of Fj 
towards a constant value. 

Equating T corr and Th y d as given above yields 



-1/2, 



1 + VI + 4ar 



(3) 
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FIG. t± '. Tp £LS £1 function of r 3 xp for different initial conditions: 
Ni = 5 • 10 4 , pi = 1.1 • 10 9 cm- 3 , T c = 50K; Ni = 4 • 10 4 , 
Pi = 3 • 10 9 cm _3 , T e = 45K; JV S = 5 ■ 10 4 , pi = 10 9 cm- 3 , 
T c = 33.3K; Ni = 8 ■ 10 4 , pi = 10 9 cm- 3 , T c = 38K; JVj = 10 5 , 
pi = 1.3 ■ 10 9 cm- 3 , T c = 33. 3K (left to right). The solid line 
is a linear fit. The error bars show the range of two to eight 
percent relative deviation between Ti and Ti for determining 



with x = T eX p/(3a) as the time when both timescales 
become equal. In fig.0J we show the time Tp when corre- 
lations start to freeze out as a function of T^ xp , where Tp~ 
is determined as the time when the relative deviation be- 
tween Tj and T{ is less than five percent for the last time. 
The linear correlation visible in the figure strongly sup- 
ports our reasoning that it is the cross-over of timescales 
that is responsible for the freeze-out of correlations. 

Thus, we may conclude that the system ultimately ap- 
proaches a non-equilibrium undercorrelated state again 
due to the correlation freeze-out described above. Still, 
the pair correlation functions can well be fitted to those of 
an equilibrium plasma at this stage (fig-EJd)), in contrast 
to the behavior at early times. This is due to the fact 
that the system went through a phase where equilibrium 
spatial correlations have developed which are preserved 
during the further evolution of the plasma. Hence, the 
system has the correlation properties of an equilibrium 
system, however "with the wrong temperature" . 

In conclusion, we have simulated an expanding ultra- 
cold neutral plasma with special attention to the forma- 
tion of ionic correlations. We have found that several 
phases can be distinguished in the evolution of the sys- 
tem. First, a quick relaxation to local equilibrium occurs, 
together with its characteristic transient oscillations of 
the ion temperature. After that, the system is close to a 
- changing — local equilibrium. In this stage, a CCP 
defined from temperature and density indeed is a mea- 
sure for correlations in the plasma. Moreover, and this 
has, to our knowledge, not been pointed out so far, the 
plasma reaches a state well inside the strongly coupled 
regime, with I 1 ; > 10. Ultimately, the timescale for equi- 
libration becomes longer than the timescale on which the 



equilibrium changes, thus the system cannot equilibrate 
anymore and correlations freeze out. Clearly, ultracold 
neutral plasmas are unique systems that evolve through 
different thermodynamical stages of non-equilibrium and 
(near)-equilibrium behavior. Their further experimental 
and theoretical study thus should provide new stimulus 
for plasma physics as well as for non-equilibrium thermo- 
dynamics. 
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